data-raw/bootstrap slopes.R

In durkeepk/sloper: Compute, extract, and visualize response-contingency slopes

library(sloper)
library(data.table)
library(tidyverse)


vars <- c("r_id", "t_id", "strength_rating", "t_measured_strength")

df <- sloper_exdat[ , vars]


rater_IDs <- "r_id"
target_IDs <- "t_id"

N_ratings <- length(unique(df[,target_IDs]))


boot_data <- data.frame(data.table(df)
                        [,.SD[sample(1:.N, N_ratings, replace = T),],
                          by = rater_IDs, ])

int_slope <- get_slopes(boot_data,
                                response = "strength_rating",
                                contingency = "t_measured_strength",
                                compress = "rater", maximal = T)

merged_data <- sloper_exdat %>%
  select(r_id, rater_formidability, rater_grip_strength, ratersex) %>%
  full_join(int_slope, by = "r_id") %>% unique()

summary(lm(slope ~ rater_formidability, data = merged_data))